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The generation of an entangled coherent state is one of the most important ingredients of quan- 
tum information processing using coherent states. Recently, numerous schemes to achieve this task 
have been proposed. In order to generate travelling-wave entangled coherent states, cross phase 
modulation, optimized by optical Kerr effect enhancement in a dense medium in an electromagnet- 
ically induced transparency (EIT) regime, seems to be very promising. In this scenario, we propose 
a fully quantized model of a double-EIT scheme recently proposed [D. Petrosyan and G. Kurizki, 
Phys. Rev. A 65, 33833 (2002)]: the quantization step is performed adopting a fully Hamiltonian 
approach. This allows us to write effective equations of motion for two interacting quantum fields of 
light that show how the dynamics of one field depends on the photon-number operator of the other. 
The preparation of a Schrodinger cat state, which is a superposition of two distinct coherent states, 
is briefly exposed. This is based on non-linear interaction via double-EIT of two light fields (initially 
prepared in coherent states) and on a detection step performed using a 50 : 50 beam splitter and 
two photodetectors. In order to show the entanglement of an entangled coherent state, we suggest 
to measure the joint quadrature variance of the field. We show that the entangled coherent states 
satisfy the sufficient condition for entanglement based on quadrature variance measurement. We 
also show how robust our scheme is against a low detection efficiency of homodyne detectors. 

PACS numbers; 42.50.Dv, 42.50.Gy, 03.67.-a, 42.65.-k 



I. INTRODUCTION 

The generation of a Schrodinger cat state [Q, which 
is a superposition of two distinct coherent states, and an 
entangled coherent state, which is an entanglement of the 
coherent states, serves the first step towards quantum in- 
formation processing using coherent states . Numerous 
schemes have been proposed in order to generate such a 
kind of coherent superposition |§, |5|. Cavity quan- 
tum electrodynamics (CQED) seems to be a promising 
environment to this task |^ and, in order to investigate 
their properties with respect to decoherence, recently a 
scheme to generate a mesoscopic version of a cat state us- 
ing trapped ions has been proposed However, most 
of the suggested schemes for quantum computation using 
coherent states is based on travelling-wave fields. Yurke 
and Stoler's suggestion to produce a travelling-wave cat 
state was far from the experimental realization because of 
an extremely low efficiency and a high absorption rate of 
nonlinear Kerr interaction while the suggestion primarily 
depends on it. 

Recently, it has been proved that the interaction of 
two travelling fields of light in an atomic medium ||^, ^ is 
able to show giant Kerr non-linearities by means of the 
so-called cross phase modulation . Measured values 
of the x*-'^^ parameter are up to six orders of magnitude 
larger than usual |pd| . This can open the way toward the 
use of this kind of non-linear processes even for the very 
low photon-number case iQ. Usually, the approach to 



such processes is restricted to a semiclassical level: the 
medium is treated quantum mechanically while the inter- 
acting fields are assumed to be classical objects. Never- 
theless, a fully quantum treatment of non-linear dynam- 
ics is relevant with respect to many aspects of quantum 
information processing. For example, huge Kerr non- 
linearities can be exploited in order to perform compu- 
tation, as said above and as described in refs. to 
perform quantum teleportation of an unknown state [l^ 
or for quantum non-demolition measurements ||l5|| . In 
all these examples, a complete quantum treatment of the 
fields involved is required. 

A full quantum analysis of the cross phase modula- 
tion problem has been explicitly performed by Lukin and 
Imamogiu in ref. ]T^ , where a rather involved atomic 
system, realized by mixing two different isotopes of the 
same alkali species, has been used. In order to suggest 
a more feasible experimental realization of the process, 
Petrosyan and Kurizki suggested a modification of the 
atomic model that allows the use of just a single species 
p^ . Their analysis, however, was again semiclassical. 

In this paper, we investigate the fully quantum- 
mechanical description of adopting a completely 
Hamiltonian approach |l^. To the best of our knowl- 
edge, this method has never been used in this context. 
We envisage in a Y2Si05 crystal doped with Pr'^+ ions a 
good candidate to physically embody the atomic model 
we discuss: the scheme of the atomic energy levels, in 
this system, seems quite appropriate to be used for our 
purposes. The Pr'^+ doped Y2Si05 has been used for 
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the experimental demonstrations of EIT Ig^l a-nd giant 
Kerr nonlinearity ||2^, Using realistic values for the 
atomic parameters relative to this solid state system we 
find that a giant rate of non-linearity is obtained in our 
fully quantum-mechanical model. We derive the relative 
equatios of motion for the involved quantum fields. This 
allows us to write the interaction Hamiltonian in a form 
that explicitly depends on the photon-number operators 
of the two quantum fields. Starting from this point, we 
show how entangled coherent states and Schrodinger cat 
states are generated when the initial states of the fields 
are two independent coherent states. 

The paper is structured as follows: in Section || we 
describe the Hamiltonian approach we have chosen and 
apply it to model cross phase modulation via clcctromag- 
netically induced transparency (EIT) In Section 
we apply this method to the atomic scheme for double- 
EIT suggested in and we derive the equations of mo- 
tion for the quantized fields. Section |^is devoted to the 
generation of entangled coherent states and Schrodinger 
cat states of light. Finally, in Section 0, we describe in 
full detail a scheme for the detection of the entanglement 
in the generated entangled coherent state. The detec- 
tion scheme is based on the total variance criterion for 
continuous variable states iial. 



II. THE HAMILTONIAN METHOD 

The standard method to describe the interaction of 
electromagnetic fields in a resonant medium is to derive 
the Bloch equations for the atomic density matrix ele- 
ments, which are solved in steady state conditions. The 
solutions are, then, inserted into the Maxwell equations 
to show the propagation of the fields. However, when 
the number of fields involved in the problem is high and 
the atomic system consists of several energy levels, this 
procedure can be quite cumbersome. 

A much simpler way to derive the field equations is 
given by a fully Hamiltonian approach |]l9| . According to 
it, the polarization of the medium can be expressed as the 
partial derivative, with respect to the electric field am- 
plitude, of the averaged free-energy density of the atomic 
medium. In other words: 



dH 

dE* 



(1) 



where H is the interaction part of the Hamiltonian, E is 
the complex amplitude of the electromagnetic field and 
P is the polarization of the medium When several 
electromagnetic fields interact with the medium, H can 
be expressed, following as: 



(2) 



where Bij are the diagonal elements of the non-linear 
third-order susceptibility, responsible for the non-linear 



terms of the refractive index at frequency ujj. Thus, the 
polarization due to the j-th electromagnetic field can be 
written as: 



' ~ \ dE* 



Ndj / dH 

~ \ dn* 



(3) 



e-^i^^t-''^'^) + c.c. 



where flj is the Rabi frequency relative to the j-th field, 
dj is the dipole matrix element of the corresponding tran- 
sition and N is the density of the atomic medium. 

Introducing this equation into the Maxwell-Bloch 
equations, we get rid of the atomic variables, obtaining a 
set of equations of motion for the Rabi frequencies that, 
in the slowly varying envelope approximation (SVEA), 
reads: 



Nd^ 




7^ . Vj. (4) 



dz'^ cdt '~ 



Changing the reference frame into — z,T = t — z/c, the 
above equation can be reduced to: 




(5) 



This approach has been used in ref. |g3[ to investigate 
the problem of resonant forward four- wave mixing based 
on EIT. It is particularly convenient if an open-system 
Hamiltonian model is used to incorporate ab initio the 
decay rates of the atomic levels and if the atoms follow 
adiabatically the fields evolution. While the first condi- 
tion can be satisfied using an effective complex Hamilto- 
nian, the second point needs more explanations. 

Solving the Bloch equations that describe the atomic 
density matrix evolution, one usually invokes the so- 
called weak coupling limit: the fields that couple an ini- 
tially prepared, collective, atomic state to other states of 
the atomic model are assumed to be very weak (usually, 
there is less than one photon per atom on average). Thus, 
the probability that, after the interaction, a state differ- 
ent from the initial one is populated is very small. This 
qualifies the initial state as a stationary state and the 
system will evolve in an adiabatic fashion, following its 
dynamics. In these conditions, the averaged Hamiltonian 
that appears in Eq. (||) can be replaced by the eigenvalue 
of H that, in the limit of vanishing weak coupling fields, 
gives the energy of the initially prepared state p3[. Thus: 



dX 



2heoc \dn* 



(6) 



where A is the above cited eigenvalue of H. 

The advantages of this approach are evident: the 
knowledge of the eigenspectrum of the single atom model 
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suffices to derive the field equation of motion directly. 
The quantization of the fields is then performed in the 
canonical way, just replacing the classical field variables 
in the effective Hamiltonian represented by the explicit 
expression of A and assigning appropriate commutation 
rules to them [Q. Starting from this effective, fully 
quantized Hamiltonian, the quantum generalization of 
the equations of motion for the fields is easily derived. 
We want to stress here that adopting this Hamiltonian 
approach we do not introduce any other approximation 
with respect to the semiclassical case: we just eliminate 
the atomic variables evolution from that of the fields 
without solving the corresponding Bloch equations. 

Here, we propose an example to illustrate the power of 
the Hamiltonian approach and to show how to get huge 
non linear effects using the interaction of a field with a 
macroscopic atomic ensemble in the EIT regime. 

We refer explicitly to ref. (the atomic model is 
sketched in Fig. ||) where, using the usual semiclassical 
approach, it has been proved that giant values of the 
third-order atomic susceptibility x^^'^ can be obtained. 
This result is a consequence of the a.c. Stark shift ex- 
perienced by the assumed metastable state |3) because 
of the dispersive coupling, induced by field £'2, between 
1 3) and |4). The Hamiltonian that describes this inter- 
action, with a canonical transformation and introduc- 
ing the decay rates of states |2) and |4), in the basis 
{|1) |2) , |3) , |4)} have the effective matrix representation: 




1> 



FIG. 1: Sketch of the energy levels of the model by Schmidt 
and Imamoglu . Fields Ei and E2 are assumed to be weak 
with respect to the strong driving field with fequency ujd- 72 
and 74 are the decay rates of states |2) and |4) respectively. 
States |1) and |3) are assumed to be metastable. In condi- 
tion of two-photon Raman resonance, the ensemble appears 
transparent to field Ei that propagates inside it with a very 
slow group velocity ^ ^. A is the detuning of the |3) ^ |4) 
transition: this dispersive coupling induces a.c. Stark shift in 
the state |3). This results in a shift of the refractive index 
curve of the medium. Because of the steepness of this curve 
inside the EIT frequency window for field Ei, the value of 
x'^' is strongly enhanced [M. 



H' 




(7) 



Assuming, as in |jT^, that |f2i| , |172| ^ , ^,74, the 
secular equation for H' results in a fourth-order polyno- 
mial whose coefficients can be expanded in power series 
of |ili,2| / l^^dl- Retaining just the first significant terms 
of these expansions, the relevant eigenvalue is found to 
be: 



A 



SI 



'(A-i74)|r!rf| 



(8) 



Taking the limit |r2i|,|f22| 0, A5/ tends to zero, 
which is the energy of the initially prepared state as 
assumed in p^ . 

Having As/, the technique described in the preceding 
section can be straightforwardly applied: deriving Xsi, 
which now represents an effective interaction Hamilto- 
nian, with respect to fi^ allows us to get an expression 
for Pi, polarization of the medium at frequency wi. Ac- 
cording to Eq. we have: 



N\d 
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'|d34|' 



(A - 274) in^r 



and then: 



iVMl2|'M34|' 

eo(A-i74)|^^<i|'ft3' 



(9) 



(10) 



This is exactly the main result obtained by Schmidt and 
Imamoglu |10 . Taking the real part of the complex x*^^^ 
we get the rate of non-linearity of this process. Note 
that, differenly from the works in rcfs. 2^, here we 
do not have any x^^^ because of the assumed perfect res- 
onance in the transition |1) <-> |2) and the zero atomic 
decay rate from state |3). Measured values of the non 
linear refractive index, for this model, are of the order 
of 10~^cm^ /Watt, resulting in an enhancement of the 
Kerr effect up to six orders of magnitude with respect to 
the best measured values for the case of cold trapped Cs 
atoms fill] . 

The main result of this section has been to show that 
the chosen Hamiltonian approach is able to reproduce 
correctly the results obtained by solving the equations of 
motion for the atomic density matrix elements. Starting 
from it, we will straightforwardly derive the full quantum 
description of a model for double-EIT. 



III. CROSS PHASE MODULATION VIA A 
DOUBLE EIT EFFECT 

We refer again to Fig. |]for the details of the following 
discussion. As explained above, in the EIT regime, the 
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field El travels in the medium with a very slow group ve- 
locity (17 m/ sec in ||ll| and Abm/ sec in |22)) while E2 has 
a very high propagation velocity. Harris and Hau proved 
that the total phase shift experienced by field Ei is 
limited by the time that the faster of the two fields spends 
inside the medium. The efficiency of the non-linear in- 
teraction is, thus, strongly affected by any velocity mis- 
match. In order to get rid of this bottleneck, strategies 
to induce EIT for both Ei and E2 (doubie-EIT regime) 
have been developed. This will maximize the interaction 
time, optimizing the efficiency of the process. While the 
scheme suggested in ref. [|l6|, even if extremely stimu- 
lating, seems to be hard to be experimentally realized, 
Petrosyan and Kurizki |l^ proposed another scheme for 
double-EIT to simplify the model. Even if it implies a 
complication of the atomic energy spectrum, it appears 
simpler under a realizable point of view. The energy 
scheme is sketched in Fig. ||: it involves a six-level atomic 
configuration and four electromagnetic fields. A mag- 
netic field splits metastable triplet {|1) , |2) , |3)} by A^, 
and the excited triplet {|4) , |5) , |6)} by A[/ ^ A^. Tran- 
sition 1 2) ^ 1 5) is assumed to be forbidden, while |2) is 
resonantly coupled to states |4) and |6) by means of the 
two very weak probes Ea and Ej,, respectively. These 
two fields couple transitions |1) ^ |5) and |3) ^ |5) with 
a detuning |A| = |Ac/ - Al|. The couplings |1) ^ |4) 
and |3) ^ |6) are realized by two classical, intense fields 
of different frequencies but equal Rabi frequencies. In 
these conditions, the system divides itself into two parts. 
For the subsystem composed of |1) , |4) , |2) , |5), EIT is 
induced for field Ea while an a.c. Stark shift effect on 
state |1) is determined by E^ to generate the required 
non-linear interaction. For the subsystem composed of 
|3) , |6) , |2) , |5), an analogous discussion can be done in- 
terchanging Ea and E^. The two subsystems are related 
via the non-resonant couplings involving |5). The double- 
EIT regime is, thus, established. 

As we have discussed above, a Hamiltonian approach 
reveals its advantages when several atomic levels are in- 
volved. In these cases, even if a Maxwell-Bloch approach 
is still possible, the procedure itself is rather uncomfort- 
able. Furthermore, the generalization to a fully quantized 
version of a non-linear process can be hard to perform 

®- . . . H . 

For the system described in Fig. H, we write the Hamil- 
tonian in the interaction picture: 

H' = hA |5> (5| + h{nd |4) (1| + Qd |6) (3| 

+ fi„|4)(2|+r!b|5)(l|+r!,,|6)(2| (n) 

+ nae~^'^*\5) (3| +C.C.}. 

To show that our method is able to mimic the results 
obtained, at the semiclassical level, by the approach cho- 
sen in ref. , we appropriately change the signs in front 
of each Rabi frequency in Eq. (|l^) and we introduce 
the excited states decay rates. In ref. ||l^, the sign in 
front of each Rabi frequency is chosen according to the 
Clebsch-Gordan coefficient of the corresponding transi- 
tion. Here, this is performed in a phenomenological way 



just to match our model with the one reported there. We 
finally get the following matrix representation of H' for 
the atomic basis {|1) , |2) , |3) , |4) , |5) , |6)}: 
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(12) 

The solution of the Schrodinger equation for a state 

|<^)=^A,(i)N) (13) 
1,1 

is obtained assuming the weak field limit ,7, A 3> 
\ila.b\ and that both 7"^ and A~^ are larger than T, 
the characteristic interaction time of the applied fields 
with the atomic medium. Under these conditions, we 
can use the SVEA for the atomic probability amplitudes 
Ai{t) [i = 1,..,6) and for the field amplitudes Ea and 
i?f,: this is equivalent to require that the amplitudes of 
the applied weak fields do not change too much during 
T. Neglecting the highly oscillating terms, in such a way 




FIG. 2: This figure sliows the atomic model used to get a 
double-EIT regime for fields Ea, with frequency uja, and E^,, 
with ujb- The fields that have frequencies ujdi, Wd2 are as- 
sumed to be classical, in the sense that their intensities is 
much greater than that of Ea and Eb- The splitting Au is as- 
sumed to be different with respect to A_t. |A| = \Au — A_l| is 
the detuning of field Ea relatively to the transition |1) <-> |5) 
and of field Eb with respect to |3) <-> |5). The excited states 
decay rates are assumed to be equal to 7, for sake of simplic- 
ity. 
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that a kind of rotating wave approximation (RWA) is 
performed, the probabiHty ampHtudes reach stationary 
values. Note that this second assumption agrees with an 
adiabatic solution of the equations of motion. 

If all the atoms in the ensemble are initially prepared 
in state |2), in the weak fields limit we can take A2{t) ~ 1, 
V< < T. Ref. 1 17 shows that the atomic polarizability of 
the medium at frequency uja is given by: 



2mo7 l^bl 
{j + iA)\nd\ 



(14) 



where ao = Nao and ctq is the resonant absorption cross 
section, generaUy defined by ctq = s'l'^- Eq. (|l|) 
shows explicitly the effect of the cross phase modulation 
induced by the interaction between the two weak (but 
classical) fields: the polarizability at frequency uja, due 
to field Ea, depends on the intensity of field £'(,. Since 
a completely analogous expression holds for the polar- 
izability ab at frequency LUb, the cross effect is evident. 
Here, we are assuming that the atomic ensemble is man- 
tained at a sufficiently low temperature to discard any 
Doppler broadening. Rigorously speaking, the thermal 
distribution of the atomic velocities has an influence on 
the value of the susceptibility of the medium, that has 
to be averaged over the velocity distribution function. If 
the temperature of the sample is kept low (orders of 10^ 
nK in |jll| ) and if we adopt a co- propagating beams con- 
figuration in order to get rid of residual Doppler shifts, 
the broadening can be made small and the average can 
be avoided [ [l7| - 

Introducing aa into the equation of motion for Ea we 
get the solution: 



EaiL,t)^Ea(^0,t-^^exp!^iJ^ aadz^ 



(15) 



where L is the interaction length of the fields. Ea ex- 
periences, thus, a phase shift due to the presence of the 
second field |6[ 0. 

We now apply the Hamiltonian approach, to show 
that the results obtained solving perturbatively both the 
Maxwell equations for the fields and the equations for 
the atomic probability amplitudes can be obtained just 
by looking for the eigen-energy of the system that, for 
\^a,b\ — *■ 0, gives the energy of state |2) (that is the ini- 
tial state of the system). 

The secular equation for the matrix H' given in 
Eq. (HI) is a six-order polynomial expression whose co- 
efficients depend on the Rabi frequencies ^a,b,d- In the 
weak field limit, we use a series expansion of fla.b/^d- 
Retaining just the first orders and performing that kind 
of RWA that allows us to neglect all the very highly os- 
cillating terms, we finally get: 



2h\naf ifibl^ \nd 



(16) 



with\n\^^\nar + \nb\ +W. 

As pointed out, this equation has to be contrasted with 
that for the eigen-energy of state |2) in absence of the 
weak probes. If these fields are absent, it is easy to ver- 
ify that the atomic model shown in Fig. |^ can be mapped 
into an effective five-level system that does not include 
state 1 2). The diagonalization of the resulting Hamilto- 
nian (obtained from Eq. (12) getting rid of the second 
row and column of the matrix) shows that the states |1) 
and 1 4) are dressed by the intense field with frequency 
Wdii while the field with frequency 0Jd2 dresses the tran- 
sition 1 6) ^ 1 3). This shows that |2) is the only state 
that, in absence of weak fields but with the strong clas- 
sical ones shined on the medium, has zero energy. Since, 
for \ila,b\ — > 0, we have A — > 0, Eq. ( |l6| ) is the right 
solution. 

Assuming once more the weak field limit, we have — 
\Qd\ and the expression for A can be approximated to: 



A 



2h\na\^ lUbf 

(z7- A)|r!rf|' 



(17) 



The partial derivative of Eq. ( |17[ ) with respect to fi* 
gives us an explicit expression for the polarization of the 
medium at frequency tOa and the equation of motion for 
r^a. The latter, finally, reads: 



dQa _ 2iNao-f\nb\ 



(7 + zA)|f]rf| 



r^a, 



(18) 



that exactly corresponds to the result semi-classically ob- 
tained in ref. [ pT| . In the same way, the partial derivative 
with respect to 17^ leads to the polarizability at frequency 
oJb and to the equation of motion for the Rabi frequency 
fif,. Because of the symmetry of the system with respect 
to the two fields, we easily recognize that their group ve- 
locities are both equal to Vg ~ Ifld] /Nao'f <^ c. This 
inequality holds using the values reported in [p^ . As ex- 
plained above, this equally slow propagation of the two 
fields inside the medium optimizes the cross phase mod- 
ulation effect. 

To have an efficient non-linear process, the rate of two- 
photon-absorption has to be negligible with respect to 
the rate of non-linearity. Since the former quantity is 
proportional to the imaginary part of the polarizability 
while the latter is proportional to his real part |l^, we 
can consider the following figure of merit for the non- 
linear interaction: 



^joa} 

^{Ua} 



A 
7 ' 



(19) 



If the experimental conditions are such that A 3> 7, 
then any absorption can be neglected and the process 
can be seen just as a mutually induced phase-shift of 
the fields. Note that this is fully consistent with the 
requirement advanced in the original theory of giant Kerr 
non-linearity by Schmidt and Imamoglu ]To| . Adopting 
the values chosen in ref. p7|, an interaction length of 
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the order of cm and an interaction time of some /isec 
lead, for two focused beams Ea and E^, to a total phase 
shift (obtained integrating 3ff{aa} over the interaction 
length) that can easily reach tt. With these orders of 
magnitude, the total two-photon-absorption probability 
is smaller than 1%. 

The quantization of the fields, now, proceeds as follows: 
we replace the complex Rabi frequencies that appear in 
Eq. (^_7|) with the positive and negative frequency com- 
ponents of the corresponding field operators (that sat- 
isfy the bosonic commutation rules [f2i,17j] cx 6iji, with 

dij the Kronecker symbol, 1 the identity operator and 
i,j = a,b), multiply the expression that is thus obtained 
by the density of the atoms in the ensemble {N) and in- 
tegrate over the interaction volume V = AL, with A the 
effective cross section of the fields. Following this recipe, 
we get an effective Hamiltonian operator that describes, 
in a completely quantum picture, the non-linear interac- 
tion of two quantum fields that propagate inside a dense 
medium in condition of double-EIT: 



H, 



eff 



2hAN 
(«7 - A) 



(20) 



For the case of pulses propagating inside the non-linear 
medium, we can follow a treatment analogous to that 
developed in p7| . Thus, adopting the narrow bandwidth 
approximation and assuming a finite range of frequencies 
involved in the superpositions that build up the pulses, 
we introduce the slowly varying positive frequency oper- 
ator f2,(z,t) = d2iEk^/^^ak{t)e-'^^^-<^^>/^ and 

the analogous for r2f,(z, t), where we explicitly introduced 
the annihilation operators a^. Here, Vq is the quantiza- 
tion volume, fc is a label for the different wavelengths 
appearing in the superposition and lu^""^ is the central 
(carrier) frequency of the pulse: the narrow bandwidth 
approximation consists in assuming that the width of 
the pulses, in the frequency domain, is smaller than 
the carrier frequency itself. When the spatial integra- 
tion is carried on, assuming that the medium length is 
longer than all the wavevelengths in the pulses, the main 
contribution to Heff is due to terms as J2k '^ki^)^k{t) 
(and the same for field Ei,) which define the total pho- 
ton number operator in the pulse Na (Ni,). In the 
case of single-mode field, as a cw laser beam, the sum 
that appears in the definition of ^a,biz,t) collapses and 

na{z,t) = d24y^^ae-*(""*-'=»^). Assuming 77 » 1 

and properly collecting all the non-operatorial quantities 
into a rate of non-linearity x, we can write Eq. ( po|) as: 



with 



Heff — hxa^ab^b 



NuJaLJb\d24\'^ \d26f 

\2h^el{ij~A)\nd\^V^ 



(21) 



(22) 



where we have assumed that the interaction volume co- 
incides with the quantization one. This is the counter- 
part, specialized to the particular atomic model we have 
adopted here, of the rate of non-linearities obtained by 
Lukin and Imamoglu in jl^ . 

Let us briefly turn to some experimental details. A 
promising candidate to embody the atomic model we 
used for the double-EIT regime is a crystal of Y2Si05 
doped with Pr^+ ions (Pr:YSO) |2^, both for an interest- 
ing similarity between the energy-level scheme described 
here and that of the transition ^i?4 ^Z?2 in this crystal 
and for the possibility we have, in a solid state system, 
to limit the effect of the Doppler broadening. This solid 
state system is notable for its relatively narrow-linewidth 
EIT that, very recently, enabled the observation of ultra- 
slow group velocity (approximatively 45 m/sec) and stor- 
age of light pulses (measured delay times greater than 65 
/isec) Ig^. For the '^H^ transition at a wave- 

length of ~ 600 nm considered in |2^, the ground state 
population lifetime is of the order of minutes. The crys- 
tal sample can be taken as long as 1 mm and the laser 
beams used in the non-linear interaction can be focused, 
by a lens, to have a typical diameter of 100 /im (Full- 
Width- at-Half-Maximum). A realistic value for the ex- 
cited states decay rate 7 range between 10 and 100 kHz 
(close to the measured values for a sample at a tempera- 
ture of 5 K) that, for detuning A ~ 1 MHz and coupling 
field Rabi frequency \^ld\ ~ 1 MHz allows us to consider 
7 <C A, \fld\, in accordance with our treatment. The 
electric dipole matrix elements for the system in exam 
are three orders of magnitude weaker than those, for the 
same range of frequencies, in alkali atoms. Taking, as 
typical values for an alkali atom d ^ 10^^^ Cm, we can 
assume |c?24| — 1^261 ^ 10^'^^ Cm. Finally, the atomic 



density is taken as iV ~ 10 cm Putting these values 



in Eq. ( p2[) we find that interaction times r in the range 
of /isec allow to reach a cross-phase shift ip — xr ~ tt 
even when the intensity of the beams involved (propor- 
tional to {fia^b)) is no more than a few photons. We will 
see in the next section how important, for the purposes 
of this paper, is the ability to get a tt cross-phase shift of 
the interacting fields. 

For the case of pulses interacting in the atomic sample, 
some experimental difficulties rise. First of all, effects of 
diffraction, focusing and defocusing on Ea and E^ are 
to be considered: they are due to the transverse inten- 
sity profile of Vid that leads to variations in the radial 
refractive index experienced by the weak probes (this ef- 
fect is known in literature as electromagnetically-induced- 
focusing (EIF) [^). EIF changes the size of the weak 
beams from point to point inside the medium, thus mod- 
ifying the interaction volume (that becomes a function of 
the position) and influencing the rate of non-linearities. 
The effect, present even in the cw regime, is less con- 
trollable in the case of pulses because of the different 
frequencies involved in the propagating packets. These 
non-linear effects, indeed, seem to depend on the sign 
of the detunings between the components of the probes 
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and the frequency of the atomic transitions that they 
guide. This impHes, for example, focusing for some com- 
ponents of a probe pulse accompanied by defocusing of 
all the other components (of the same pulse) that have 
a detuning of opposite sign. The result is that different 
harmonics, in a pulse, are subject to different radial evo- 
lutions, thus complicating the control on the dynamics of 
the field itself. 

Another relevant point to be treated is that of the dif- 
ferent phase shifts acquired by the different parts of the 
cross-interacting pulses. As it can be seen solving the 
Heisenberg equations of motion for the field operators 
Ea and Eh interacting as described in Eq. (^o|), and as it 
is proved in , the evolution of the probe fields is given 
by: 



EaAL.t) = ^a..(0,i')e{'*<»^°'*')^''-(°'*'^> (23) 



with t' = t — L jvg and x the rate of non-linearity specific 
for this case (that reduces to Eq. (22) when we consider 
cw fields). Note that this solution is the exact quantum 
analogous of what has been found, in the semiclassical 
approach, in Eq. (p^. Thus, the phase shift experienced 
by pulse a (b) depends on the total number of photons in 
pulse b (a) at the earlier time t' and this number changes 
with the amplitude of Et (Ea). The effect is that different 
parts of a propagating pulse acquire different phase shifts 
with respect to each other. It modifies, in essence, the rel- 
ative phase relations between the pulse components and 
distorts the pulse shape. Spectral width enlargement, 
for example, can be a detrimental consequence for the 
purely dispersive propagation inside the atomic medium: 
if some harmonics of the evolving pulses exit from the 
EIT transparency window, they will be strongly absorbed 
by the no-more-transparent medium. Usually, a way to 
bypass this kind of problem is to arrange the pulses to be 
within the EIT window since their entrance into the non- 
linear medium, properly choosing their shape and spec- 
tral width. For adiabatic evolution and for an optically 
thick medium, the pulses will last in the non absorptive 
region of the refractive index of the medium through all 
the interaction time |^ . But this result has been proved 
just for a simple atomic A system and it needs a deeper 
analysis when we refer to the atomic energy scheme we 
describe in this paper. In our opinion, all these points 
have to be further investigated, even seeking for an ex- 
perimental verification of their real influence, in order 
to have a complete comprehension of the kind of con- 
trol we can reach for the cross-phase interaction of two 
weak fields in a medium that exhibits giant Kerr non- 
linearities. 



For the sake of simplicity, but without lacking the nec- 
essary experimental realism, in what follows we treat just 
the case of single cw mode propagation and, from now on, 
we proceed taking in consideration the effective Hamil- 
tonian in Eq. (pT|). 



IV. SCHRODINGER CAT STATES 
GENERATION 



In this section we want to apply the results obtained 
with our fully quantized picture of the cross phase modu- 
lations induced by double-EIT in order to show that this 
specific system can be used to produce a Schrodinger cat 
state of a single mode of field. 

Given the interaction Hamiltonian (|2^), we derive the 
equations of motion for the annihilation operators a and 
b. Adopting the usual notation ha = a^a, hb — wb for 
the photon-number operators of the two fields, these are: 



dta = -^x[a,?^.a^^^,] = -ixafib 
dtb = -ix[b,hahb] = -ix^na 



(24) 



which, arranging the time scale in such a way that t = 
is the instant in which the interaction starts, lead to the 
time-dependent operators: 



a-out{t) 



a(0) bout{t) = e-'^'''-b{0). (25) 



Note that, in the approximation of a pure cross shift pro- 
cess, ha,b do not evolve because of the non-absorption 
character of the interaction considered. Defining ip = 
the above evolution is attributed to the action of the 
unitary time-evolution operator U {(p) = e"*'''"""'' on the 
field operators. We want to specialize the present analy- 
sis to the case in which the initial state of the two inter- 
acting fields is \'ip{0))ab — l"^)a® \^)b' where the coherent 
state \a)^ — Da{a)\0)a with the displacement operator 
Da{ct) — exp(aa^ — Q*a) The coherent state |7)f, 

has been defined likewise. 

The evolution of the initial state by means of U {ip) can 
be expressed as follows: 



\m)ab=Da{a)Db{l)\0)a<»\0)b, 



(26) 



where we have defined the time-dependent displacement 
operators: 



Da{a) = U{cp)Da{a)U-\ip) = eC""^-"*'^) 
^,(7) = U{^)Dbh)U-\v) - e(^^'-^*^) 



(27) 



with at = p-'vnafn^-lgtvAam and ^ ^-zvA^fnl^f^tvAaA,^ 
Explicitly using the operator expansion theorem |]3l|| , we 
obtain: 



that is: 




_ (jb^e-'^"--j'be"^<'-) 



(28) 



(29) 
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Introducing Eq. into Eq. (| 

dependent state of the two fields: 



|), we get the time- 



Da{a)D,{j)\0)^(g>\0), 



|0),®|7)6, (30) 



where the Campbell-Baker-Haussdorff theorem and 
the fact that a coherent state is an eigenstate of the an- 
nihilation operator have been used. Using the represen- 
tation of coherent states in the Fock number states we 
have: 



n,0 



\ae 



-tip' 



/ a 



(31) 



If we are experimentally able to set the interaction time 
and the value of the rate of non-linearity so that ip = t: 
then, splitting the sum in Eq. into one over the odd 
values of n and one over the even ones, the evolved state 
of the two initially non-interacting fields takes the form 

W^/X))ab « l«)a {l7) + 1-7)};, + |-«)a {l7) - |-7)}b • 

(32) 

This is a particular expression for an entangled coher- 
ent state: it can be reduced to the more familiar form 
1*^)0 l7)b + \~'^)a l~7)& unitarily acting on the subsystem 
b. To prove the entanglement, we have to show the corre- 
lation of the fields of modes a and b as we unitarily trans- 
form, gradually, from [7)^ to |7)j, ± |— 7)^. However, this 
involves another non-linear interaction. We thus discuss 
an indirect way to prove the production of the entangled 
coherent state. 

In Eq. (^2|), linear superpositions of the coherent states 
|7)j and |— 7)^ are Schrodinger cat states: 



2i 



|2j), 



l7)h- |-7)& oc 



v2j + l 



(33) 



|2j + 1) 



which are sometimes called as the even and odd coherent 
states. As shown in Eq. (|3^), a time-controlled inter- 
action of two fields in initially prepared coherent states 
results in an entangled state: if we properly normalize 
state |V'(^/x))a6 look to the reduced density op- 

erator of the field Ea alone (tracing over b mode) , we find 
an incoherent mixture of two coherent states which are 
out of phase by tt . 

Suppose we are, somehow, able to discern the state in 
which field Ea is, after its crossing through our model 
of highly non-linear medium. Because of the entangled 
structure of Eq. (^, the state of field Eb is projected 
onto one of the equally weighted superpositions of 7r-out- 
of-phase coherent states [7)^ ± |— 7)^: it collapses into a 
Schrodinger cat state. The point, now, is to show how to 
reliably discern the state of field Ea- 



To achieve this target, we need a 50 : 50 beam splitter 
(BS). After passing through a beam splitter, two coherent 
input fields |a)|/3) become [p3[: 



Ba 



a + P 



V2 



-a + (3 



V2 



(34) 



where Bac = e^'^"'^ with a (a) and c (c) as the 
input (output) modes of the beam splitter. 

Thus, taking P = a and referring to Fig. ^ we have the 
following read-out scheme: if the input mode a is in the 
state \a)a, then Detector 1 will click, revealing that some 
photons arrived at it while Detector 2 will not click. In 
this case, the field mode b will be projected in the even 
coherent state. In the opposite occasion, the field mode 
b will be in the odd coherent state. Of course, there is a 
possibility to have both the detectors not to click. In this 
case, we do not know where the mode a is so we have to 
repeat the experiment till we have one detector to click. 



field mode a 



Detector 1 
field mode a 




Detector 2 



FIG. 3: Detection scheme to infer the state of the field mode 
a. In this figure it is shown the symbol used to represent 
a photodetector and that for the 50:50 beam splitter (BS). 
According to Eq. (p4|). Detector 1 (2) clicks just if the state 
of mode a is |q:)^ Jj^a)^). Using this read-out scheme, we 
have the possibility to generate an even or an odd coherent 
state of the field mode b. 



V. DETECTION OF A CAT STATE 

The quantum state of an electromagnetic field is com- 
pletely and unambigously determined once all the mo- 
ments of its statistical distribution function are deter- 
mined To get this task, we need the density matrix 
p of the state of the field: it contains all the available 
information about the given quantum state and, having 
/3, the most complete statistical description of the sys- 
tem can be performed. On the other hand, it is possible 
to demonstrate a biunique correspondence between the 
density matrix p and the Wigner function of a radiation 
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mode so that, determining the latter, we are able to 
fully identify the quantum state of the field. This cor- 
respondence is useful under an observable point of view: 
numerous theoretical schemes for the reconstruction of 
the Wigner function of a radiation mode have been pro- 
posed. They are essentially based on heterodyne or ho- 
modyne detection p5[ . The mathematical manipulation 
of the data collected by the homodyne detector allows 
also for the reconstruction of p ^ . Recently, this ho- 
modyne tomographic technique has been experimentally 
realized p8| . This approach, however, may be limited by 
imperfect detection efficiency. 

The procedure described above could certainly be used 
to determine if our scheme for the non-linear interaction 
of two weak fields inside a Kerr-enhanced medium has 
been able to generate a Schrodinger cat state. However, 
as its complete and satisfactory treatment is far beyond 
our purposes, we want to suggest alternative ways to 
show the signature of the quantum nature of the state 
generated. Recently, a sufhcient criterion for the insepa- 
rability of a continuous variable state has been proposed 
by Duan et al. They proved that, for a separable 

continuous variable state, the total variance of a couple 
of suitably defined conjugate operators respects a lower 
bound imposed by the uncertainty relation. For an en- 
tangled state, on the other hand, this bound can be vio- 
lated, providing a sufficient criterion for inseparability of 
the state. More precisely: assume a continuous variable 
state composed of subsystems 1 and 2 and take the cou- 
ple of operators, for both subsystems, "Pij [i — 1, 2) 



such that 



A' V- 



— i6ji. Because of this commutation 



rule, the simultaneous measurement of both the opera- 
tors is affected by the Heisenberg uncertainty principle. 
Following fl^ , we take an arbitrary, non zero, real pa- 
rameter q and define the pair of collective operators: 



u=\q\Xi + -X2 

q 

V = \q\Vi - -V2. 



(35) 



For any separable state, it is possible to show that the 
following inequality holds Km : 



S^{{Mlf) + {{Mf)>q^ + q- 



(36) 



with <^(Au)^) and <^(At))^) the average variances, calcu- 
lated over the state of the system, of u and v respectively. 
In what follows, we will often refer to the total variance S 
as to the separability function. It is worthwhile to stress 
that the criterion is just a sufhcient condition for insepa- 
rability. The established bound could, thus, be exceeded 
even by a continuous variable state that is not separable. 
We will give an explicit example of this possibility in the 
following discussion. 

The criterion by Duan et al. leaves a certain freedom 
in the choice of the couple of operators to use we need to 
construct u and v. The only limitation imposed is that 
they have to be conjugate operators. 



In problems related to the investigation of the quantum 
state of light, it is often useful to consider the phase-space 
of the quadratures: 



(37) 



So defined, the quadratures of a field are a couple of con- 
jugate operators, since [x,p\ = i. Operatively, their prob- 
ability distribution can be reconstructed by a homodyne 
detection scheme ||^, ^: the input mode to measure is 
mixed, at a 50 : 50 beam splitter, with the coherent state 
of a local oscillator whose phase is directly controllable 
during the experiment and whose intensity is so high to 
consider it classically. The difference between the num- 
ber of photons in the output modes is, then, measured 
by two photodetectors (one for each output mode). The 
data collected are proportional to the expectation value 
of the operator: 



(38) 



that, for i} = coincides with x, while for -d = tt/2 
is equal to p |3^. The couple of quadratures {xa,Pa}, 
{xb,Pb} relative to the field modes in Eq. ( ^2| ) are thus 
well suitable to construct operators u and v for the joint 
system of modes a and b. What we need, to calculate the 
total variance function 5, is just the variance of single 
quadrature (such as ((Aia.b)^) for example) and the ex- 
pectation value of correlations as (xaXb) or (paPb), mea- 
surable collecting the coincidences at the detectors Q . 

We take q = 1, so that the bound value for the sep- 
arability of the input state is 2. Unfortunately, the cal- 
culation of the separability function S for the entangled 
superposition in Eq. (p2) leads to S* > 2, whatever are the 
amplitudes a and 7. As we stressed above, this certainly 
does not mean that the state is separable but just makes 
this criterion unsuitable to reveal the entangled nature 
of the investigated state. Because the behaviour of S is 
very state-dependent, to bypass the negative result we 
obtain using Eq. (|3|) directly, we propose the following 
strategy. Assume we generated, as described at the end 
of section [V, an even coherent state of mode b. We then 
mix it with the vacuum of mode c at a 50 : 50 BS. Spe- 
cializing the general rule of a beam splitter in Eq. ( |3^ ) 
for /? = 0, the joint state of the BS output modes can be 
written as: 



with TV = 



7 

V2 



7 
'V2 



7 

'V2 



7 



, ^ This is an entangled coherent 

V2(l+e-2^') 

state of two field modes: note that, using this procedure, 
this output state can be obtained only if the radiation 
mode b is in an even coherent state. More generally, to 
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describe the state of mode 6, we can take the density 
matrix: 

Pb = ^{|7> (7I + |-7> (-7I + c(|7) (-7I + |-7> (7l)}b , 

(40) 

with < |c| < 1 and A a normahzation constant. For c = 
0, the state is a statistical mixture; c = 1 (c = — 1) gives 
us the density matrix of an even (odd) coherent state 
while, the general case in which < c < 1 corresponds 
to a non optimal generation of the Schrodinger cat state. 
After the action of the 50 : 50 BS, the two output modes 
h and c (whose initial density matrix was |0)^ (0|) will be 
described by: 



7 -7 \ 1-1 
V2' 72/ \V2' V2 



-7 7 \ / -7 7 
^/2' 72/ \V2' V2 



+c 



-7 7 
V2' n/2 



7 -7 
V2' V2 



7 -7 
%/2' x/2 



z2 

V2' V2 



be 

_ (41) 
It is evident that state (|39| ) can be retrieved from this 
density matrix only if c = 1 (we have the case of an odd 
coherent state in input when c = —1). 

The quadrature statistics of the two radiation modes 
are measured by two homodyne detectors, as described 
above. The collected data allow to estimate the exper- 
imental total variance for the input state. This has to 
be contrasted with what we get calculating the separa- 
bility function. We thus test, by means of the sufficient 
criterion for inseparability, the entanglement properties 
of state (^9|) in order to infer the state of the radiation 
mode h. Even if, by this means, we are certainly not sure 
that the state of modes h and c is exactly the entangled 
coherent state in Eq. (|3^), we will show the entangled 
nature of the detected state. 



If we calculate S for the case of state we get: 



5, 



perfect 



2U 



2e-2^ 



1 



-272 



(42) 



A plot of Sperject as a function of the amplitude 7 is 
given in Fig. H: the total variance function stays below 
the bound S = 2 just until 7 ~ 2 (we have S = 1.995, 
for 7 = 2) and except 7 = 0. This gives the signature of 
the entangled nature of the low intensity state we have 
generated. The dip of the separability function shows a 
maximum deviation from the 5 = 2 bound equal to 28% 
for 7 ~ 0.8 and is still about 5% for 7 = 1.6. 

Thus, the criterion seems to work sufficiently well for 
low values of the amplitude of the input coherent states. 
This is because reducing the value of the amplitude, the 
entangled coherent state becomes similar to a Gaussian 
state, that is a class of states for which the total vari- 
ance criterion for inseparability works very efficiently. In 
Fig. (|^), the separability function ( p2| ) is compared to 
what is obtained taking c = and c = 0.5 (that is mid- 
way between a statistical mixture and a perfectly gener- 
ated entangled coherent state). The largest effect of an 
imperfect generation of the input is to reduce the dip of 



the non-separability well until, for the statistical mixture 
corresponding to c = 0, the criterion fails for whatever 
value of the input amplitude. The information we obtain 



Sperfec 




FIG. 4: Plot of Sperfect = ((Au)^) + {(A.vf) as a function 
of the coherent state amplitude 7. The sufficient criterion 
for inseparability states that, for any separable continuous 
variable state, this total variance function has to be bound 
from below by a value that, for the calculations shown in 
this paper, is equal to 2. The effect of different values of the 
parameter c in the density matrix p'^ _ is studied: the dot- 
dashed curve is for c = 0, corresponding to to the case of a 
statistical mixture. The dashed curve is for c = 0.5 while the 
solid curve represents the case of a perfectly generated even 
coherent state of field mode b. In this latter case, for the state 
considered in Eq. (^), the bound is beaten just until 7 ~ 2. 

with the detection scheme we are describing give us a suf- 
ficient insight into the entangled nature of the state (|39|). 

We now show the robustness of the detection scheme 
for entanglement with respect to the unavoidable homo- 
dyne detector inefficiencies. In what follows we attribute 
them essentially to losses. Other possible sources of error 
are the dark-counts at the photodetcctors, whose entity 
depend on the width of the detection window. Their ef- 
fects can be made negligible taking the detection time 
short with respect to the characteristic time of a dark- 
count occurrence. Finally, we neglect the noise intro- 
duced by the local oscillators supposing their amplitudes 
to be well stabilized. 

The effect of loss on the collected data can be mod- 
eled replacing a real, inefficient, homodyne detector with 
a beam splitter (transmittivity 77) followed by a perfect 
homodyne detector In Fig. |5|the experimental appa- 
ratus that considers inefficient homodyne detectors ac- 
cording to the above equivalent model is sketched. Each 
beam splitter BS,,, that mixes a mode of the signal to 
measure with a vacuum state, transmits the signal with 
probability i] and reflects it with probability 1 — rj. The 
amount of reflected input field gives a measure of the 
losses of the detection apparatus, so that the quantum 
efficiency of this latter can be identified with rj itself. To 
avoid stray light to enter in the perfectly efficient homo- 
dyne detector, we can immagine that the reflected output 
field enters a blackbody absorber. Mathematically, this is 
equivalent to tracing out the degrees of freedom relative 
to the absorbed field mode. What the perfect homodyne 
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detectors measures, thus, is the statistic of the quadra- 
tures of the transmitted beams, irrespective of the state 
of the reflected ones. 




FIG. 5: Sketch of the proposed scheme to infer the value 
of the total variance for the entangled state of modes b and 
c when the inefficiency of homodyne detectors is taken into 
account. With Di and D2 we indicate perfect homodyne de- 
tectors that have to measure the quadratures of modes d and 
e. These are the fields transmitted (with probability -q) by 
two beam splitters, BS,,, whose effect models the inefficiency 
of the homodyne detectors. The beams reflected by BS,, are 
absorbed by blackbodies. The local oscillators are prepared 
in very intense coherent states. This allows to treat them 
as classical fields. Their phases are experimentally controlled 
and adjustable. BSi and BSi are 50 : 50 beam splitters. 



Simperfect 




y 



0.5 i 1.5 2 2.5 3 



FIG. 6: Behavior of the separability function for state (p9[), 
as a function of the amplitude 7, when imperfection in the 
homodyne detection are taken into account. Here, the dot- 
dashed curve is for detection efficiency rj = 0, the dashed one 
is for ?7 = 0.4 and the solid curve is for rj = 0.8. For lower 77, 
the minimum values of Simperfect shift toward higher 7 values: 
this is because the lower is the efficiency of the homodyne 
detectors, the more state (^) resembles a Gaussian state. 

The calculation of the total variance for the quadra- 
tures of modes b, c when the detectors have an equal 



Probability Distribution 
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FIG. 7: Plot of P^.-q{x), for an input field in the even coherent 
state with 7 = 1.5. The effect of different values of rj is 
investigated: the dot-dashed curve is for 77 = 0.1, the dashed 
one is for 77 = 0.4, the solid one is for t] = 0.8. 

quantum efficiency 77 leads to what is shown in Fig. 
The separability function seems to keep its functional 
features even in the case of imperfect detection. How- 
ever, some differences, with respect to Sperfect, appear. 
As the quantum efficiency of the homodyne detectors re- 
duces, the mimimum of the total variance shifts toward 
higher values of 7 (see Fig. ^), enlarging the range of 
amplitudes for which the separability function has values 
below 2. 

This effect can be explained directly referring to the 
probability that the homodyne detector measures a value 
X of the quadrature x. Without loss of generality, but 
avoiding the lengthy calculation relative to state (^), we 
investigate an even coherent state. We let this state to be 
mixed with a vacuum field at a BS of transmittivity 77. As 
before, the resulting transmitted mode is then measured 
by a perfect homodyne detector while the reflected field 
is traced out. We finally get the following probability 
distribution for the in-phase quadrature x: 



(43) 

It shows that P-y^ji.oi^) is composed of two Gaussian 
hills centered at zLy/2rjj plus a third term that, for 7 
sufficiently large, is negligible with respect to the oth- 
ers (even for 7 = 2 it is negligible). The larger is the 
amplitude 7, the larger the Gaussian hills are apart 
However, if an efficiency 77 < 100% is considered, the dis- 
tance between the Gaussian hills for low values of 7 is 
strongly reduced and the joint state appears more and 
more similar to a Gaussian state. To observe well sepa- 
rated Gaussian bells we have to go toward larger 7. It is 
possible to show that the same effect is observed in the 
probability distribution of an entangled coherent state as 
Eq. (|39|), with the suitable modifications due to the dif- 
ferent nature of the state. This explains the behaviour 
of the separability function for an imperfect homodyne 
detector. 

However, the success of the total variance criterion for 
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low values of the quantum homodyne inefficiency has not 
to be considered positively. The imperfections in the de- 
tection device tend to hide the features of the input state 
that appears as a Gaussian state even if, in reality, it is 
well far to be Gaussian. This means that, in order to 
reliably test the entangled nature of the input state, a 
sufficiently large ry has to be considered. Realistic values 
of the homodyne detector efficiency range between 0.6 
and 0.85. In particular, according to our calculations, 
r] = 0.8 seems to be a good trade off between the relia- 
bility of the total variance criterion and the realism of an 
experimentally achievable quantum homodyne efficiency. 

From the above analysis, it appears that for 7 larger 
than 7 ~ 2, it is not possible to get any information 
about the entanglement in the input state of modes b 
and c from the total variance of quadrature amplitudes. 
We have to look for other detection strategies. 

If we still refer to the case of a generated even coher- 
ent state in Eq. (|3^), the characteristic of the superpo- 
sition of 1 7) J and |— 7)}, is well revealed by homodyne 
measurements, as pointed out by Yurke and Stoler 
As commented above, the amplitude of the components 
of an even coherent state can be inferred directly from 
the probability distribution of the in-phase quadrature. 
If instead of x we measure the statistical distribution for 
■d = 7r/2 in Eq. (38), the following function has to be 
found: 



Probability Distribution 



1 



-2(l-r,)7- 



cos(2-y/2777a;)| 



-272 



(44) 

This probability oscillates and is moduled, at the same 
time, by a Gaussian function (solid curve in Fig. ^). The 
oscillations, with a frequency dependent on 7 as well as 
on ?7, are an evidence of the quantum interferences be- 
tween the two out-of-phase components of the coherent 
superposition in the even coherent state. It is, indeed, 
straightforward to prove that, if instead of a linear com- 
bination of 17)1 and |— 7)t, the input state was the mixed 
state of Eq. (^0|) with c — 0, then this oscillatory behav- 
ior would be absent. Thus, a homodyne detection scheme 
is able to discern a Schrodinger cat state from a mixed 
state of two coherent states which are mutually n out of 
phase Hj. 

It is important to stress the effect, on the oscillatory 
pattern, due to the Gaussian modulation: as it is evident 
from Fig. ^, even a small reduction of rj leads, for a small 
amplitude coherent state, to the disappearance of the 
oscillations. Mathematically, this is due to the factor 
1 — r] that governs the width of the Gaussian function 
in front of the oscillatory term: the smaller 77 is, the 
more rapidly the Gaussian function goes to zero, washing 
out the interference pattern. The effect is, clearly, more 
evident for large values of 7. This means that this scheme 
is very vulnerable to the detection inefficiencies. 

An alternative approach to the investigation of the co- 
herences in the quantum superposition of [7)^ and 1—7)^, 
when the total variance criterion fails, can be the fol- 




FIG. 8: Probability distribution for a detected homodyne cur- 
rent with 1} — ■k/2. In this plot, 7 = 2 while is scanned: the 
dashed curve is for 77 = 0.4, the dot-dashed one is for 77 = 0.9 
while the solide curve is for r] = 0.95. 



lowing. We prep are, once more, the entangled coherent 
state in Eq. (|39|) using an even coherent state as input. 
Similarly to the technique used by Kwiat et al. in pl| , 
we operate unitarily on one of the components of \4')i^- 
Then, we detect the coincidences between the counts col- 
lected by two photodetectors facing the output modes c 
and g (see Fig. ^). 

In the restricted Hilbert space spanned by the orthogo- 
nal even and odd coherent states, the rotation of a generic 



state 1,5)^ = ^1^)- + 



B 



"^/~ ''^^ performed by 
a displacement operator Di{iO) {6 E M) that, acting on 
I (5)^, transforms it into: 



|A) = Ae^ 



V2 



+ ^^ 



Be VI 



If we take 9 <^ l/V^, on the Bloch-sphere of the re- 
stricted Hilbert space that we are treating, the above 
equation approximates well the state R^ica) \S)i, with 

R^{a) the rotation operator, by an angle a/2 — \/2^9, 
around the Bloch-sphere z-axis. Experimentally, such a 
rotation can be accomplished superimposing mode b on 
the coherent state \i£) of an auxiliary mode /, using a 
beam splitter with a high-transmission coefficient T. As 
shown in , the transformation operated by the T ^ 1 
beam splitter (BS2 in Fig. ||) gives an entangled state of 
two output field modes. The reduced density matrix that 
describes the state of the output mode g only can then 
be approximated by that of a displaced state according 
to pg = b~^{iE\/^~T)pibl{iEVT^) where pi is the 
density matrix of the input field |^ . 

If we vary £ in such a way that 9 = £\/l — T <C 7/\/2, 
that is: 



£ 



2^2(1 -r)7' 



(46) 



then, just controlling the amplitude of this ancillary 
mode we can appropriately rotate the state of mode b 
by a/2. 
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FIG. 9: The figure sketches the apparatus to detect the gener- 
ated state of mode b. As before, the input even (odd) coherent 
state and the vacuum state of the auxiliary field mode c enter 
the 50 : 50 beam splitter BSi and give rise to an entangled 
coherent state of the output modes b and c. While mode c is 
directly detected by the photodetector A, mode b is rotated 
by the device shown in the dashed box. It consists of the high 
transmittivity {T 1) BS2 that superimposes a coherent 
state \i£) to the field mode b. The transformation thus real- 
ized approximates well that of a displacement operator that 
rotates the field in b (see Eq. ([l7|)). The trasformed mode is, 
then, sent to detector B. BothTand B are two On/Off pho- 
todetectors that are able to discriminate the vacuum from any 
incoming photon. They are not able, however, to reveal the 
number of incident photons. The coincidences of the counts 
at the two detectors are measured. 



Following this lines, the entangled coherent state in 
Eq. (^9|), after the interaction of b with / at the rotation 
beam splitter BS2, becomes: 



7 

V2 



-7 
V2 



-7 
V2 



7 

V2 



(47) 



The output modes c and g are then sent to detector A 
and B respectively, where the photo-count coincidences 
are revealed. We model the photodetectors as Geiger-like 
On/Off Photo-Detectors {0/0 PD). These particular de- 
vices just discriminate the vacuum from an input with 
any photons, irrespective of what the photon number is 
p3[ . Again, losses make the quantum efficiency of the an 
0/0 PD not optimal, so that some of the coincidences 
will be missed. Quantum mechanically, an inefficient 
0/0 PD can be described by means of a suitably defined 
Positive- Operator- Valued- Measure (POVM), that is a set 
of diagonalizable, having positive-eigenvalues, projection 
operators. The appropriate choice for this case is the 



click 



{i = A, B) (48) 



cLick 



We assume the same efficiency 77 for both the detectors. 
In the above definitions 1 — ?y is the probability that a 
single detector lacks to reveal a photon and the sum over 
all the photon populations takes in account the impossi- 
bility to distinguish the photon number in the incident 
field. Within the POVM formalism, the probability that 
the 0/0 PD click within the same detection window is 
given by: 



P, 



W.(7,^,a) = Trg,[\^)^,mn^^l^{r^)n^^l^{r^)] 




(49) 

with C ~ ^e"^ N"^ . Let us look to the asymptotic be- 
haviour of the above probability: for 77 ^ 0, that is for 
highly inefficient detectors, we get P2ciicks{l. O1 a) ^ 0. 
If the limit 77 1 is instead considered, Eq. (|49| ) shows 

that P2ciicksil, 1, a) ^ 2 ^1 -f e^'^^ cosa^ . The oscillat- 
ing part in P2ciicks{'^ , l,a) is the evidence of the coher- 
ences established in the even coherent state. Unfortu- 
nately, it has an exponential pre-factor that can cause 
experimental difficulties (for 7 ~ 2 and 77 ~ 0.75, a vis- 
ibility of about 1% is achieved). Furthermore, possible 
errors connected to the rotation of a Schrodinger cat state 
have to be considered. 

Once the entangled nature of the input state |<&)j,g 
is recognized, however, the inference about the state of 
field mode 6, generated by the non-linear interaction via 
double-EIT, proceeds along the lines outlined before. 

From the above discussion, it is clear that the main 
problems related to the described schemes rely on the 
quantum efficiency of the detection devices. As we have 
seen, imperfections in the homodyne detectors, for exam- 
ple, are responsible for the disappearance of the oscilla- 
tion pattern in the probability distribution relative to the 
out-of-phase quadrature. For the case of a Schrodinger 
cat state, this results in a loss of information relative to 
the quantum superposition of two coherent states. It is 
worthwhile to stress that for high value of rj, the oscilla- 
tions in Pj.r1.-K/2ix) survive even for large values of 7. 

The hard task represented by the detection of an even 
(odd) coherent state or, more generally, of an entangled 
coherent state can be made easier by the recently im- 
proved sensitivity of the available detectors. Values of 
?7 ^ 0.8 are realistic for a detector operating in a Geiger- 
mode and very sensitive photo-detectors, with an effi- 
ciency as high as ~ 0.93, have been developed. They 
however require a low temperature to properly oper- 
ate El. 
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VI. CONCLUSIONS 

We have proposed a fully quantized picture of the 
model for double-EIT recently discussed by Petrosyan 
and Kurizki |l^ . Our approach is based on a full Hamil- 
tonian method. This allows to bypass the analytical so- 
lution of the equations of motion of the atomic density 
matrix elements. It simplifies the computational prob- 
lems related to systems that involve many atomic energy 
levels coupled by electromagnetic fields. 

As shown for the first time in [p^ , the non- linear in- 
teraction of two beams of light that pass through a dense 
atomic medium is optimized by a double-EIT regime in 
which both the fields propagate with a strongly reduced, 
equal, group velocity. The atomic model proposed in ||l^ 
seems to be a good candidate to the experimental realiza- 
tion of such a physical condition. In our fully quantum 
version of this model, we have shown that the quantum 
dynamics of one of the interacting beams is dramatically 
dependent on the intensity of the second and viceversa. 
Our results are consistent with those reported in . We 
suggested the PriYSO crystal as a possible candidate 
to physically embody the discussed atomic model. 

Starting from these results, we have written an effec- 
tive interaction Hamiltonian that, when the field is ini- 
tially prepared in two coherent states, leads to the pos- 
sibility to generate entangled coherent states and even 
or odd coherent states (Schrodinger cat states). In or- 
der to investigate the quantum features of a generated 
even (odd) coherent state, we described a scheme able to 



measure the total variance function of the involved modes 
quadratures ^ . For low amplitudes of the investigated 
entangled coherent state, the value of this total variance 
is well below the bound imposed by the inseparability 
criterion suggested in ||l^. We also investigated the ef- 
fects of an imperfect non-linear process and of detection 
inefficiencies on the above results. We found our scheme 
robust against homodyne detection losses. For the cases 
in which the total variance criterion is unuseful, we sug- 
gested a detection scheme based on the unitary manipu- 
lation of one of the modes in the entangled coherent state 
followed by the detection of coincidences of photo-counts. 
Wc discussed possible sources of errors and experimental 
difficulties we expect to appear in a real experiment, find- 
ing that an optimal detection protocol basically requires 
high efficiency of the detectors. 
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